On 

o\ 

-<— > 

o 
O 

(N 

S 1 
Pk 

CZ2 

Renormalization Group Studies of Vertex Models 

' m ' 

>v 

Saibal Mitra 
Supervisor: Prof. Dr. H. van Beijeren 
<N i ITF-UU-99/04 

m ; 9th February 2008 

O 

o\ 

On 

czj 

O 

• i— i 

> 



Contents 



1 Introduction 



1 . 1 Summary 



3 

4 



Definition of the staggered F-model 5 

2.1 The six- vertex model 5 

2.2 Six-vertex models and SOS models] 6 

2.3 Roughening transition in the F- model 7 



3 Renormalization group equations for sine-Gordon type models 8 



3.1 Effective Hamiltonians for the staggered F-model 

3.2 The renormalization group transformation .... 



3.3 Renormalization 10 



3.4 Cumulant expansion 11 



3.5 Diagrammatic expansion 12 



3.6 Evaluation of diagrams 16 



Applications of the renormalization group equations! 20 

4.1 The case of the staggered F-model 21 

... 22 

. . 25 



4.2 



4.3 



Singular part of the free energy of the staggered F-model 
The case % < | In (2 - y/2)\ 



5 Expansions about free fermion models 



5.3 The Pfaffian method 29 



28 

5.1 Definition of free fermion models 28 

5.2 Baxter's solution of the staggered F-model 28 



5.4 Calculation of the free energy 31 

5.5 Singular part of the free energy 33 



5.6 Perturbation theory! 35 

36 



5.7 First order computation for the staggered F- model 



5.8 



Singular behavior in the vicinity of the free fermion line 39 

40 



5.9 Linked cluster expansion 



A Notations and conventions 



47 

A.l Vectors] 47 

A. 2 Summation convention 47 



A. 3 Multi-indices 47 



A. 4 Derivatives 48 



A. 5 Distributions 48 



B 



Gaussian correlations 49 

eTi Correlation function for hS 1 ' 1 49 

B.2 The height-height correlation function 50 



C Calculation of 



-TT7T 



51 



Chapter 1 

Introduction 



In this thesis we use renormalization group methods to study the critical be- 
haviour of the staggered F-model. The staggered F-model, defined in chapter 
pi can be used as a model of a facet of a BCC crystal in the (100) direction. 
We first use known exact results to map the staggered F-model to a sine- 



Gordon type model (defined in section 3.1 ), and study the renormalization group 



equations for this model using momentum shell integration techniques. The map 
to the sine-Gordon model is constructed using an exact result on the long range 
part of the height-height correlation function of the F-model (i.e. the staggered 
F-model at zero staggered field) pO| . The results we obtain are the phase 
diagram of the staggered F-model and the leading singularity in the free energy. 
To get more results, e.g. the next to leading singularity in the free energy, we 
need a larger source of information than is available in the form of the long range 
part of the height-height correlation function. It turns out that the free fcrmion 
method, upon which Baxter's original solution is based, is flexible enough to 
admit a perturbative expansion about the free fermion line. By calculating the 
singular part of the free energy perturbatively about the free fermion line, one 
can construct a map from the staggered F-model to the sine-Gordon model 
by demanding that it correctly reproduces the singular part of the free energy. 
The idea thus is to use the mapping to the sine-Gordon type model as an 
extrapolation technique. To make this approach practical we: 

1. Develop in chapter |3| a simple diagrammatic method to find the renormal- 
ization group equations for a given sine-Gordon type model. This method 
is based on a combination of functional Feynman rules and the operator 
product expansion. 



2. Rewrite in section 5.9 the perturbative expansion about the free fermion 



line as a linked cluster expansion. 

Although these two results make it possible to construct a map to a sine-Gordon 
type model in a systematic way, the actual construction of this map is beyond 
the scope of this thesis. We do, however, explicitly calculate the first order 
correction to Baxter's result. This allows us to verify results which previously 
could only be obtained using renormalization group arguments. 



1.1 Summary 

In chapter we introduce the staggered F-model and discuss the equivalence 
with the BCSOS model. 

In chapter y a systematic method is derived to generate the renormalization 
group equations 

In chapter |j we first discuss the application of the renormalization group in 
the calculation of critical exponents. We then proceed to obtain the phase 
diagram of the staggered F-model and calculate the leading singularity 
in the free energy by using the information present in the form of the 
asymptotic form of the height-height correlation function 

In chapter g we first present a derivation of Baxter's exact solution. We then 
perturbatively lift the free fermion condition. This allows one to write the 
free energy of the staggered F-model as a perturbative expansion about 
the free fermion line. Baxter's exact solution can thus be seen as the zeroth 
order term in this expansion. We explicitely calculate the first order term. 
To facilitate the computation of the higher order terms we derive a linked 



cluster method in section 5.9 



Chapter 2 

Definition of the staggered 
F-model 



In this chapter we shall introduce the six-vertex model, of which the staggered 
F-model is a special case. Some known results are discussed. 



2.1 The six- vertex model 

The six-vertex model can be defined as follows: place arrows on the bonds of 
a square lattice so that there are two arrows pointing into each vertex. Six 
types of vertices can arise (hence the name of the model). These vertices are 



shown in fig. 2.1. By giving each vertex- type a (position-dependent) energy 



the model is defined. These models were first introduced to study ferroelectric 
systems. Later it was shown that six- vertex models can be mapped to solid-on- 
solid (SOS) models M. Only a few of these models can be solved exactly. These 
include the free fermion models |q,|39| and models that can can be solved using 
a (generalized) Bethe Ansatz ||,g,|2l[ |23j. To define the staggered F-model, we 
divide the lattice into two sublattices A and B, such that the nearest neighbor 
of an A vertex is a B vertex. The vertex energies are chosen as indicated in 
fig. 2.1. When the the staggered field (s) vanishes the model reduces to the 



F-model, which has been solved by Lie b p2|| . For nonzero staggered field the 
model can be solved when j3e = \ In (2) j|. 



5 
±s 



G 



Figure 2.1: The six vertices and their energies. The upper and lower signs correspond 
to the two sublattices. 
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Figure 2.2: An arrow configuration together with the corresponding height function. 



2.2 Six-vertex models and SOS models 



We now proceed to show how six-vertex models are related to SOS models. 
First we introduce a dual lattice. Each bond of the dual lattice now crosses an 
arrow placed on one of the bonds of the original lattice. By rotating this arrow 
90° clockwise and placing it on the corresponding bond of the dual lattice, we 
obtain an arrow configuration on the dual lattice. A height function (h) is now 
defined by demanding that h (x) — h (y) + 1 if an arrow points from y to x. 
By fixing the height at one particular point, the height at each point of the 
dual lattice is defined unambiguously. See j3| for more details. The fact that 
the height difference between nearest neighbors is always ±1 makes six- vertex 
models ideal models for crystal surfaces of BCC crystals in the (100) direction. 
The class of SOS models to which six-vertex models are mapped is also known 
as body centered solid on solid models (BCSOS models). In fig. 2/2 



an arrow 



configuration on a lattice together with the corresponding height function on 
the dual lattice is shown. 



2.3 Roughening transition in the F-model 

According to p2[ a phase transition of Kosterlitz-Thouless type takes place in 
the F-model at inverse temperature /3e = In (2). If (3e > In (2) the crystal surface 
as described by the F-model is smooth. In this case the height-height correla- 
tion function G (r) = / (h (r) — h (0)) ) decays exponentially with increasing r. 



When one takes /3e < In (2), the surface is in a rough phase. It can be shown 
that Q 

G W = TT^ ^T ln M ( 2 - X ) 

7T arccos (1 — | exp (2pe)J 

The logarithmic divergence of the correlation function at large distances is 
caused by thermal fluctuations in the local height of the surface with arbitrary 
long wavelengths. Note that for e > the F-model has a twofold degenerate 
ground state consisting of vertex 5 on one sublattice and vertex 6 on the other 
sublattice. By introducing a staggered field this degeneracy is lifted. It has 
been shown [E7J that in a nonzero staggered field the F-model is in a smooth 
phase for positive e. 



Chapter 3 

Renormalization group 
equations for sine-Gordon 
type models 



In this chapter we will introduce the sine-Gordon type Hamiltonian and then 
show how renormalization group equations can be obtained for such models. 
First a cut-off procedure will be introduced to define the theory. Renormaliza- 
tion is carried out by first integrating over some of the degrees of freedom of 
the model. The model, when formulated in terms of the remaining degrees of 
freedom, will look like the original model with a lower cut-off. Finally a scale 
transformation will restore the original cut-off. 

3.1 Effective Hamiltonians for the staggered F- 
model 

Since the staggered F-model can be interpreted as a solid-on-solid model (see 
section [2^2| ), it is natural to introduce a field h, that describes the height of a 
surface. The Hamiltonian density of this field must possess the same symmetries 
as the staggered F-model. In particular we must have: 

F(h + l,0s) = F(h,-/3s) (3.1) 

F(h) = F(-h) (3.2) 

Here s is the staggered field, and we have assumed that the ground state of 
the staggered F-model (for (3e > and (3s ^ 0) corresponds to ft. = in the 



sine- Gordon model. From (3.1) it follows that 



F(h + 2,(3s)=F(h,(3s) (3.3) 



(p.2j) and (3.3) lead us to the Hamiltonian density: 



F (h, dih, dijh, ... , (3e, (3s) — 2_] D n (dih, dijh, ... , j3e, (3s) cos (mrh) 

n=0 (3.4) 



Here D n is an unknown function of its arguments. According to ( p.l[ ) we have 

D n (dAd v h,... ,/3e,-/3s) = (-l) n D n (d i h,d ij h,... ,/3e,/3s) 

(3.5) 

3.2 The renormalization group transformation 

We will rewrite the Hamiltonian fl3.4|) as 



00 f d 2 x 
H = \J / —j- exp(irnrh)D n (adih, a 2 dijh, a 3 dijkh 7 . . . ) (3-6) 

—00 

We can think of the constant a as the "lattice constant" of the original micro- 
scopic Hamiltonian. In this original model \ would be the density of degrees 
of freedom. The effective Hamiltonian <pM should have the same density of de- 
grees of freedom. The constant a appears in the Hamiltonian as a consequence 
of replacing summations by integrals and finite differences by partial derivatives. 
We will define the Fourier transform of the field h(x) as 

h (k) = —j= j d 2 xh{x) exp (— ik ■ x) (3-7) 

v V J 



Here V is the volume of the system. h(x) can then be written as 

h(x) — —^= y^ h (k) exp (ik ■ x) (3.8) 

We now define a cut-off by introducing a set (S) of allowed A:- values. We assume 
that the set S has the property: 

keS^-keS (3.9) 

The density of A:- values is written as , ^ P (fc). The function P (fe) will be 
called a cut-off function. We shall assume that the cut-off is chosen such that 
P (0) = 1 and all derivatives of P (k) are zero at k = 0. If the volume V is 
chosen large enough, a summation over S can be replaced by an integral: 



E 

fees 



F(k) = V f-^P(k)F(k) 

J (Z7r) 



(3.10) 



provided that the function F does not correlate with the characteristic function 
of S. In case such correlations do exist we have to replace P (k) by the char- 
acteristic function of the set S, which we denote as P c (k). In general we thus 
have 

Y,F(k) = V f^P c (k)F(k) (3.11) 

fees J ( 27r ) 

The value of a now follows by requiring \ to be the number of degrees of 
freedom per unit volume: 

1 f d 2 k , N , 

-=■ = / ~ p £0 (3.12) 

a 2 J (2tt) 2 



We will denote the set of all allowed functions by S. S is the set of all hnite linear 
combinations of the functions e lk ' x with k 6 S. Note that we have h (k) = if 
he S and k (jL S. 

We now define the partition function as: 

Z= f Dhe H (3.13) 



Where the measure Dh on S is defined as 

dh(k) 
kes 



Dh ^Y[^±R(k) (3.14) 



The function R (k) which occurs in the definition of the measure has to be chosen 
such that the free energy of the exactly soluble Gaussian model is consistent with 
the renormalization group equation for the free energy. Although the correct 
choice of R (k) is important for a consistent description of the theory, it turns 
out that its effect is equivalent to adding a constant term independent of any 
couplings to the Hamiltonian and hence doesn't influence the dependence of the 
free energy on the couplings. 

3.3 Renormalization 

The renormalized Hamiltonian is obtained from (|ij) by using the Wil son-K ogut 
momentum shell integration technique J2§, f37|. We will integrate ( |3.13| ) over 
some of the degrees of freedom, leaving us with an effective Hamiltonian (H) 
with a lower cut-off. Next a scale transformation will restore the original cut-off 
and yield the renormalized Hamiltonian (Hr). 

We must now specify precisely the degrees of freedom we have to integrate 
over. Since the renormalized Hamiltonian (Hr) has the same cut-off function 
P (k) as the original Hamiltonian (H), and since it is obtained from the effective 
Hamiltonian (H) after a scale transformation, H has to have a cut-off function 
of the form P (Ik). In terms of I the scale transformation becomes x — > l~ l x. We 
thus have to construct a set S^ of allowed fc-values for H, such that £W q g 
and S^ corresponds to the cut-off function P (Ik). The complement of S**- 1 ^ in 
S, denoted as S^ 2 ' , contains the degrees of freedom we have to integrate over. 
We thus have to split the set S of k- values into two disjoint sets S^ and S^ 2 \ 
This can be done as follows: We decide to put the points k € S and —keS 
in S^ with probability %&. S^ is defined as S^ = S- S^ . The cut-off 

function for S^ will be denoted as P^ , is thus given by 

F (2) = P (k) - P (Ik) (3.15) 

We now construct the spaces S^ and S^ analogous to S: S^ is defined as 
the set of all finite linear combinations of the functions e lk ' x with k G S^K We 
now have 

£ = s«es (2) (3.i6) 

The projection of a h € S on S^ and S^ will be denoted by h^ 1 ' respectively 
h,( 2 \ The first step in the Wilson-Kogut renormalization scheme is to integrate 



10 



over the field h^ 2 \ After this integration one obtains an effective Hamiltonian 
H which depends on h^\ The final step is to restore the original cut-off by a 
length rescaling: 

x/ = l~ 1 x (3.17) 

The renormalized field h R is defined as: 

h R (x/) = h^(x) (3.18) 

and the renormalized Hamiltonian Hr is defined as: 

H R (h R ) = H (h™) (3.19) 



3.4 Cumulant expansion 

The integration over the field h^ 2 ' is performed after an expansion about the 



Gaussian model. We rewrite our Hamiltonian (3.6) as 



H = H g +X (3.20) 

where H g is a Gaussian interaction and X is a perturbation. H g may be split 
into a Gaussian interaction for hP^ and h^ 2 \ denoted as H^- 1 ' respectively H^ 



9 ~ 


-ij{Vhfd 2 x = - 


-*£*=* * a IM*)l 2 


= 


-§E fce s<i)fc 2 IMfc)l 2 


~ 2 £fces< 2 > ^ I' 1 Wl 


= 


-lj(Vh^) 2 d 2 x 


-lf(VhW) 2 d 2 x 



H^+HW (3.21) 



Note that for a given Hamiltonian the representation (3.20) is not unique be- 
cause one may choose to include a Gaussian term in the perturbation X as well. 
Such a freedom of choice can sometimes be exploited in first order calculations 
to improve the accuracy of calculations (see Q ) . 
We define the measure DhP^ by 

[ Dh (2) F(h) s Les^ DhF ( h ) = IU k es^dh(k)F(h) 

J Les^ DheH<2 '' JU k es^dh(k)e^(H( 2 )(h( 2 ))l 322) 

where F{h) is an arbitrary function of h. The Gaussian average of a function 
F over the field h^ can now be written as 



(F (h)) = / Dh (2) F (h) exp (h {2) (h {2) )) (3.23) 

We now define the effective Hamiltonian H W- 1 '") as follows: 

exp (h (h {1) )) =K I DhW exp (H) (3.24) 



Here if is a constant. To determine Hr one simply has to rescale H (see (3.17), 
(p . 1 8 ) and (|3.19[)). To fix the constant K, one has to express Hr and H in the 



11 



same functional form and then require the constant terms to be equal. From 
flOol ), ( |3T2l| ), (|J23|) and ( gjj ) it follows 



i? = In (X) + iJ (1) + In (exp (X)> (3.25) 

To second order in X, ( 3.25J ) can be written as 



H = In {K) + tf « + (X) + i ((X - (X)) 2 ) + . . . (3.26) 

This expansion is known as the cumulant expansion. For the general form of 
this expansion, see fl4||. 



3.5 Diagrammatic expansion 

It turns out that the terms in the cumulant expansion can be represented as 
amplitudes of Feynman-diagrams. In these diagrams the correlation function of 



the field h^ plays the role of the propagator. In section B.l we show that it 
takes the form: 

exp (ik ■ x) 1 f dk ( 2 ) it ex P ( J ^ ' x ) 



i v exp {ik ■ x) _ i r d\ (2) 



fe6 s(=) n J J ^"' a (3.27) 

where P c is the characteristic function of the set S^ . The amplitudes of 
Feynman-diagrams we will encounter later can be expressed as integrals of prod- 
ucts of propagators. We have to be carefull with replacing Pc by P^ 1 in such 
cases. E.g. we have 

d 2 k pW(fc) 



f ,2 r„, s-,2 1 f d 2 k P [Z >(k) , nnn . 

j^{^ x) f = - 2 j-— 2 ^l (3.28) 



It is not difficult to see that G (0) is universal: 

G(0) = ^-ln(Z) (3.29) 

We now consider the case of an infinitesimal cut-off change: 

l~ 1 = l-t (3.30) 

( |j.l7 ) becomes 

x' = (l-e)x (3.31) 

We now associate e with an infinitesimal increase in a rescaling parameter t. The 
renormalization process then generates one parameter families of Hamiltonians 
H (t) . The renormalization group equations can then be written as 

— — = coefficient of e in Hr (3.32) 

The parameter t is related to a length transformation: 

x (t) = e- l x (0) (3.33) 
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Instead of the Hamiltonian it is often more convenient to write the renormal- 
ization group equations in terms of the Hamiltonian density. We shall denote 
the effective Hamiltonian density corresponding to the effective HamiltonianiJ 
asF: 

H (>)) = I d 2 xF (hP , fi,fcW) (3.34) 

The renormalizcd Hamiltonian density, denoted as F, can thus be expressed in 
terms of F by rewriting ( |3.34| ) in terms of h R : 

H R (h R ) = H (/i«) = J d 2 xF {h^\dih^) 

= Jd 2 x'(l + 2e)P(h R ,(l-e)d l h R ---) .^ 

where in the last line we used the transformation x' = (1 — e) X and h R (x 1 ) = 
h^ (x). The renormalized Hamiltonian density (Fr) can thus be expressed as 

Fr (h R ) = (1 + 2e) F (h R , (1 - e) d % h R • • • ) (3.36) 

The renormalization group equations can thus be expressed as 

dF 

— = coefficient of e in Fr (3.37) 

at 

We now proceed with the derivation of the Fcynman-rules for the cumulant 



expansion ( 3.26 ). It is convenient to derive these rules first for the term (e x ). 
We shall see that In (e x ) is obtained by summing over connected diagrams only. 
Let F (h, dih, dijh, . . . ) be the non-Gaussian part of the Hamiltonian density in 
([3.6]). We can then write: 

fc=i \fc=i / (3.38) 



We can evaluate ( 3.38 ) by writing F, considered as a function of the field h and 
its derivatives, as a Fourier integral. We will define a Fourier transform of F as 
follows: 

F \li In In ' j ■ ■ ■ ) = J 57 J 1 li ~2$T J 1 Uj ~^T ■ ■■ C\ qq\ 

F (h, dih, d^h, ...) c-»b*+7.a,fc+7«a«fc...] V'™> 

The integrals in ( |3.39| ) are from — oo to oo. F can now be written as 

F (h, dih, d l3 h, ...)=/ d-y J ]Ji d ^ I ILj <*Ty ■■■ , > 

F ( 7 , 7l , jij,...) ^hh+^h+^o^h...] ^- 4U J 



The next step is to substitute the representation ( 3.40]) for the Hamiltonian 



density in (3.38). To facilitate this, it is convenient to introduce multi-indices. 



The term in the exponent in (3.40) can be written as 

oo 

7^ + !>!,.. aA,...,^ (3-41) 



fc=i 
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A tuple of k indices, as in the summation in ( |3.41 ), can be treated as a single 
index. Such an index is called a multi-index. A tuple of k indices will be written 



as (k). We can thus rewrite (3.41) as 

oo 



(3.42) 



fc=0 



Note that repeated multi-indices are only summed over while ke eping the num- 
ber of indices contained in the multi-index constant. See section A. 3 for all the 



conventions on multi-indices. Inserting ( 3.40 ) in ( 3.38 ) gives 
1 



(X n ) 



-i/ (nM7«}j (n/CM)) (e^^£> 



Note that the term in the exponent in ( 0.43 ) can be interpreted as the action of a 
distribution (i.e. a linear functional) on the field h. The action of a distribution 
T on a function h is denoted as Th. See section A. 5 of the appendix for a precise 
definition of distributions. The problem is thus to evaluate 

(e tTh ) (3.44) 

for a general distribution T . In chapter O of the appendix it is shown that 



iTh\ _ p -%T x T v G{x-y) iThW 



(e lTh ) 



(3.45) 



Here T denotes a distribution, T x and T y act asT on G (x — y) considered as a 
function of x respectively y (x and y are thus "dummy" -variables). 

We now have to expand (3.45) in powers of the propagator, and substitute 
the result in (3.43). The distribution T in ( p. 45 ) is defined as follows: First we 
define the distribution T"' as 



fe=0 

The distribution T is then defined as 

n 

j 1 — y^ j>Cj) 



(3.46) 



(3.47) 



3=1 



The Li order term in the propagator in the integrand of (|3.43|) becomes 



1 (~lf 
n! 2 L L\ 



[T x T y G(x-y)} L lf[F({^}) 



iTh 



(i) 



(3.48) 



\3 = l 



Each term in the expansion of [T x T y G (x — y)] can be represented diagram- 
matically. We first perform a trivial step: 



[T x T y G (x - y)] L = H T x T y G (x - y) 



(3.49) 
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Each term in the expansion of the product can be represented diagrammatically 
as follows. Draw the N points Xj . If we choose from the p term in the product 
the term 'y, \d, -J from TL and the term ^,{8, I from T„, we draw an oriented 

! {m) (m) x '(n) (n) y' 

line from x r to a: s , we label the line with the value of p, and at the points x r 
and x s we put the labels (m) respectively (n) on the line. We repeat this for 
all values of p. There is now a one to one correspondence between the set of all 
possible terms in the expansion of the product and the set of labelled diagrams. 
The amplitude of a labelled diagram is obtained by inserti ng th e appropriate 



product of the 7's and the derivatives of the propagators in ( p.48| ). We see that 
the integrals over the 7's result in a factor 



1 d r F 



* r d(d {mi} h) ■■■d(d( mr )ti) 



(3.50) 

h=hW 



for each vertex where r lines, labelled by (mi) . . . (m r ), come together. The 
pro duct of the factors -7 at each vertex will precisely cancel the factor (—1) 
in ( fr.48 ), because each propagator gives rise to two factors - and there are L 



propagators. We can simplify matters further by omitting all labels, except the 
multi-indices at both ends of each propagator, in a labelled diagram. The am- 
plitude of such a Feynman diagram is given by the sum of all the corresponding 
labelled diagrams. It is convenient to define a propagator Gt n \< m ) (p — q) as 

G (n)i(m) (p - g) = dl:[ x dl% y G (x - y) (3.51) 

Since all labelled diagrams corresponding to a nonvanishing Feynman diagram 
make identical contributions, we simply have to multiply the amplitude of one 
diagram by the number of ways of labeling a Feynman diagram, to obtain its 
amplitude (relabeling the vertices will change the amplitude of a diagram, but 
when integrated over all positions of the vertices, all diagrams obtained from 
each other by a relabeling of the vertices will, of course, make identical contribu- 
tions) . This amplitude then has to be integrated over all the Xj . We shall denote 
the number of ways of orienting the propagators, labeling the propagators and 
the vertices by respectively N±, N2 and N3. Since the multi- indices have to be 
summed over, two labelings of the propagators will not be considered distinct 
if the only difference is a permutation of the multi-indices. Two labelings of 
the vertices are considered distinct if it is not possible to transform one labeling 
into the other by a relabeling of the propagators. We then have 

Ni = 2 L - k (3.52) 

where k is the number of lines that have both there ends connected to the same 
vertex, 

where the product is over all ordered pairs of vertices, and k r denotes the number 
of propagators connecting the pair r and 

n' 

N 3 = - (3.54) 

15 



where S is the order of the symmetry group of the Feynman diagram. Using 
( pl8| ), ( pO] ), ( |j5j] ), ( pT52| ), (|H|) and (|j34|), we see that the Feynman rules 



for (e x ) are as follows 



1. To compute the contribution that is n order in X and Lr order in 
i, draw all topological distinct Feynman diagrams with n vertices and L 
lines. 

2. Label both ends of each line by arbitrary multi- indices. 

3. For each vertex there is a term: 

d r F 



d {d( mi )h) ■ ■ ■ d (d {mr) h) 



(3.55) 
h=h<v 



where the (to.;) are the multi- indices on the lines at the vertex and the 
derivative is evaluated at the coordinates of the vertex. 

4. Each line labelled with the multi-indices (m) and (n) corresponds to the 
propagator G(„),( m ) (p - q): 

G {n)Am) (p-q) = dff hx d$ )y G (x - y) (3.56) 

where p and q are the coordinates of the vertices connected by the line. 

5. For each line that has both its ends connected to the same point there is 
a factor \ . 

6. For each pair of vertices connected by k lines there is a factor ^y. 

7. There is a factor 4, where S is the order of the symmetry group acting 
on the vertices of the diagram. 

8. Integrate over all coordinates of the vertices, and sum over all multi- 
indices. 

We will now show that ln(e x ) is precisely the sum of all connected diagrams. 
We assume that all connected diagrams are enumerated in some arbitrary order. 

tin 

Let Ci be the amplitude of the i connected diagram. Using the above Feynman 
rules, we can write: 

<expW) = En^ = nE^exp E^ (3.57) 



3.6 Evaluation of diagrams 



There are two diagrams, see fig. 3.1, contributing to the first order cumulant. 
Using the Feynman rules derived above it is a simple matter to evaluate the 
amplitude of these diagrams. In the case of a Hamiltonian density F (h, dih • • • ) 
one obtains to first order the effective Hamiltonian density F 

d 2 F 



f(^),^( 1 )...)=f(^),^w...) + 1e 



d( k )hd {l) h 



G (l),(m) (°) 

h=hW (3.58) 



1G 




Figure 3.1: The two Feynman diagrams corresponding to the first order cumulant. 



where the sum over I and m is from to oo. To obtain the renormalization 
group equation for the Hamiltonian density from this, we have to perform a 



rescaling x — > (1 — e)x (see (3.36)) and use (3.37). These equations yield the 



first order renormalization group equation for the Hamiltonian density: 



§— f>«^+^ a * F G " Mo 



h d ( d ^ h ) (fc) ^d(d (k) h)d(d {k) h) e (359) 



The sum over k and I is again from to oo. The quantity G(h\n\ (0) is universal 
(i.e. independent of the form of the cut-off function P (k)) when k = I = or 
k + I = 2. It is not difficult to derive the result: 

G mm) (0) = (-1)^ ±- ^ Al+m C mm) (3.60) 

Here A n is zero if n is odd, else we have 

A n = d|fc||fcr 1 (P(|fc|)-P(|(l + e )fc|)) (3.61) 

Jo 

In particular we have: 

A = e 
h 



A ? 47r • 



(3.62) 



The tensor Cn\t m \ is a contraction operator. For an arbitrary tensor Tn\ t m \, 
Tn\i m \Cn\ i m ) is the sum of all contractions of the indices contained in (I) and 
(to) We can write Cn\ r m \ explicitly as a sum of products of Kronecker delta's: 

C {l)>(m) = 22 S *(ii)Mi2) ' ' ' ^(ii +m _i),7r(i i+m ) (3-63) 



are thus i+m . terms in the sum 



The sum is over all noncquivalent permutations of the indices i\ ■ ■ ■ ii+ m - There 

(l+m)\ 

We now proceed with the evaluation of the higher order cumulants. Ac- 
cording to the Feynman rules the n order contribution to the renormalizcd 
Hamiltonian is an rt-fold integral over the volume of a product of propagators 
and functions of the field h^> . We want to replace such an expression by a sin- 
gle integral over the volume, thus obtaining a contribution to the Hamiltonian 
density. We write the amplitude A of a diagram as 

1 = I d 2 x 1 ---d 2 x n P{x 1 ---x n )D(h ( - 1 ' ) {x 1 )---h w {x n )) 

V ' (3.64) 
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Eigenoperator 


Scale dimension 


1 





cos (nh) 


7T 

4.7 


a 2 (Vh) z - i 


2 


(a 2 (V/i) 2 - i ) cos (tt/i) 


— 4-2 


cos (2irh) 


7T 
■7 



Table 3.1: Some eigenoperators and their scale dimensions relative to the Gaussian 
Hamiltonian H g = — ^ J d 2 x (V/i) . 



Here P (#i • • • x„) is the product of propagators and D (hy> (x\) ■ ■ ■ hS 1 ' (x n )) 
denotes the product of derivatives of Hamiltonian densities. It is now tempt- 
ing to perform n — 1 of the n integrations in (3.64) by Taylor-expanding the 
field about one of the points x\ . . .x n (it doesn't matter which integrations are 
performed because different choices are related by a partial integration). The 
problem with this approach is that it assumes that the field M 1 ' is analytic. 
In reality one should expect a Taylor-expansion of the field to converge only 
in a region the size of a (3.12), because a is the distance between independent 
degrees of freedom of the field. 

A better way to proceed is to use the so-called operator product expansion. 
Before we explain how this works we will first introduce some new terminology. 
A local operator is a term in the Hamiltonian that depends only on the field in 
one point. The Hamiltonian density evaluated at a certain point is an example 
of a local operator. The first order renormalization group equations for local 
operators is almost identical to that of the Hamiltonian density. If O (h (x)) is 
a local operator then 



dO 

dt 



E* 



dO 



i d(d ik) h) 



%^+oE 



d 2 



G 



(k),(0 



k.l 



9 (d(k)h) d (d( k )h) 



(3.65) 



The only difference with ( 3.59| ) is that the term 20 doesn't appear on the r.h.s. 
of this equation. Eigenoperators are local operators that renormalize as 



dO 

~dt 



-XO 



(3.66) 



A is called the scale dimension of the operator O. In table 3.1 we have listed a 



few eigenoperators with their scale dimensions. By solving ( 3.66 ) for all eigen- 
operators one obtains a complete set of operators. All eigenoperators can be 
written as a multinomial of derivatives of the field h multiplied by cos (nirh) 
with n an integer. We shall call an operator even (odd) if n is even (odd). We 
can now expand any local operator in this set of eigenoperators. A product of 
operators localized at different points can be considered to be local if the points 
lie close to each other. This product can then be expanded in eigenoperators 
localized at one point. It is clear that this expansion, known as the opera tor 
product expansion, can be used to replace D {hS 1 ^ (x\) ■ ■ ■ h^ (x n )) in ( 3.64 ) by 

a sum of operators localized at the point x\. Suppose all eigenoperators are 

th 
enumerated by an index n. The scale dimension of the n operator will be 
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denoted as A„. We can thus put 



D (h (xi) •■■h (x n )) = ^2c k (x!--- x n ) O k (h (xi)) (3.67) 

k 



Note that we have replaced the field hf- 1 ' by h. To apply (3.67) to (3.64) a 
rescaling must thus be performed first. It is important to note that (J3.67 ) 



is an identity in the sense that the Hamiltonian to which the l.h.s. is added 
may be identified with the Hamiltonian to which the r.h.s. is added. Since 
D (h (xi) • • • h (x n )) is a sum of products of eigenoperators located at the points 
Xi • • -x n , all we need to know are the functions Ci t ...i nt k (xi • • ■ x n ) (operator 
product expansion coefficients ) in the expansion 

0»i (h (xi)) • ■ • O in (h (x„)) = ^ c»i ■•■»„,* ( x i ■ ■ ■ x n) O k (h (xi)) 

k (3.68) 

The operator product expansion coefficients can be determined as follows. One 
dema nds that the replacement of the product of eigenoperators according to 
([3.68) commutes with a renormalization. One then obtains an equation relating 



Cii—i n ,k ( x i ' ' 'X n ) to Ci 1 ...i n .k \jy • • • i)- with I the rescaling factor involved in 
the renormalization. Now, when one takes Xi = x% = • • • = X n the operator 
product expansion is trivial. The functions Ci^...i nt k { x i ' ' ' x n) can thus be de- 
termined by taking the limit I — ► oo. Note that that since the renormalization 
has to be carried out perturbatively one obtains the operator product expan- 
sion coefficients as an expansion in the non-Gaussian couplings. It is thus very 
straightforward to find the operator product expansion coefficients to zeroth or- 
der. Higher order contributions to the operator product expansion coefficients 
will again involve nontrivial integrations over Feynman-diagrams. These dia- 
grams must again be evaluated using the operator product expansion. E.g. to 
find the renormalization group equations to second order one has to deal with 



expressions as in (3.64) with n — 2. Since the function D is already of second 
order one only has to work out the operator product expansion to zeroth order, 
which is straightforward. To third order one has to calculate the operator prod- 



uct expansion coefficients in ( 3.68| ) with n = 3 to zeroth order and the operator 



product expansion coefficients with n = 2 to first order. The latter ones involve 
Feynman-diagrams in which the two operators are connected to one of the other 
operators in the Hamiltonian. These diagrams can be evaluated by again using 



the operator product expansion (3.68), but now with n — 3 and only to zeroth 
order. It is clear that repeated use of the operator product expansion allows 
one to obtain the renormalization group equations to any order. 

in the next chapter, we are going to apply the theory to find the phase 
diagram of the staggered F-model. 
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Chapter 4 

Applications of the 
renormalization group 
equations 



Once the renormalization group equations are known it is a simple matter to 
obtain the singular part of the free energy. In this section we shall first derive 
the renormalization group equation for the free energy and then proceed to show 
how the singular part of the free energy is obtained from it. 
It is convenient to rewrite the Hamiltonian as 

d 2 x {Vhf + £>„ / ~4o n (h (x)) (4.1) 



Here the O n are the eigenoperators defined by equations ( 3.65 ) and ( 3.66 ). We 
write the renormalization group equations as 

j °° 

-± = (2 - A,) Vt + J2 Yl X i,ii-ikVii ■■■Vik ( 4 - 2 ) 

k— 2 i\---ik 

Note that the \i are the scale dimensions of the operators. Also note that 
the Gaussian coupling is kept constant under renormalization. This is possible 
because a 2 (Vft.) — 4 is an eigenoperator. To find out how the singularity in the 
free energy is related to the A's in this equation, we must first find the relation 
between the free energy of the original and the renormalized system. Note that 
the renormalized Hamiltonian satisfies the relation 

exp™ = / ^P-Dh (4.3) 

where Zr and Z are the partition functions for the renormalized respectively 
the original system. Combining (3.22) with (J3.24 ) and (4.3) we get: 



§ = f 
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Here Z g is the partition function of the Gaussian model: 



z a = 



/ exp (h {2) ) Dh (4.5) 



Let U be the constant contribution to Hr from In (e xp (X )). Since the total 
constant contribution to Hr is zero, it follows from ( |3.25 ) that In (if) = — U. 
( |4.4| ) can then be written 

F=(l-2e)F R + ^ + ±ln(Z g ) (4.6) 

where F and Fr are the free energy densities times —j3 for the original respec- 
tively renormalized system. 

dF 

— -2F + c + c' = (4.7) 

Here c is the coefficient of e in y and d is the coefficient of e in ^ In (Z g ). Since 
d only depends on j, which is kept constant under renormalization, the effect of 
this term is to shift the free energy by a constant amount. We are thus allowed 
to ignore this term. 



4.1 The case of the staggered F-model 

When j3s — and f3e < In (2), it is known that the F-model renormalizcs to the 
Gaussian model 1l3, Ilf3, Il7|l . For the latter we have 



H = 2 / d x{ y h > ( 4 - 8 ) 

The Gaussian coupling j is a known function of the temperature of the F-model: 

J = \ arccos (l - i e 2 ^ (4.9) 

( |4.S| ) is valid when /3e < In (2) and is obtained as follows: The long range 

part of the height-height correlation function R(r) = ( (h (r) — h (0)) ) of both 

models show a logarithmic behaviour, and is thus invariant under horizontal 
scaling. This means that the amplitude of the height-height correlation function 
is invariant under a renormalization. For the F-model one finds |J40| 

R M (T^-i ^7mY ln M ( 4 - 10 ) 

7r arccos (1 — ^ exp (2pe)J 

In case of the Gaussian model one finds (see section |B.2| ) 

R(r) ln(r) (4.11) 

nj 

Equatin g th e amplitudes of both correlation functions then leads to the identi- 
fication Q. 
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We expect that when (3s ~ 0, we may replace (3.4) by a Hamiltonian of the 
form ( |4.1| ). Then because of (BJj) the y n multiplying even (odd) operators will 
be even (odd) functions of (3 s. We now assume that the y n in ( f4.l| ) are analytic 
in some neighborhood of (3 s — 0. This implies that the y n corresponding to odd 
operators are O (/3s). Of all operators, the operator 0\ (h) = cos (nh) has the 
lowest scale dimension: 

Al = I. (4.12) 

Between 3 — \ and 3 = \ this is the only relevant operator (i.e. an initially 
infinitesimal yi increases exponentially under renormalization) . Below j = ? 
there are no relevant operators and above j = § the operator cos (2irh) also 
becomes relevant. Because the coupling y\ becomes proportional to (3s in the 
limit (3s — ► 0, we expect the staggered F-model with an infinitesimal staggered 
field to be in a different phase than at zero staggered field for those values for (3e 



that correspond to a value for j between 3 — \ and 3 — \- According to (4.9) 
this is for (3e in the interval i In (2 — \/2) < (3e < In (2). Note that at the lower 
boundary (3e is negative: \ In (2 — y/2) « —0.2674 At zero staggered field the 
logarithmic behaviour of the height-height correlation function indicates that 
the surface is in a rough phase. If the staggered field is turned on the model no 
longer renormalizes to a Gaussian model. If the staggered field is chosen small 
enough we expect that under a renormalization the model will renormalize first 
to a model of the form 



// = / d 2 x -± (Vhf + V - 1 
I a 



-^(Vhy + ^cos(nh) 



(4.13) 



with a small value for y\, but as we renormalize further the coupling y\ will 
increase. Since the effect of the operator cos (irh) in the Hamiltonian is to favour 
even values of h, we expect to be in a smooth phas e. B elow (3e — | In (2 — \/2) 



we still expect that the model will renormalize to (3T) but as we renormalize 
further y\ will renormalize to zero. We are then left with a purely Gaussian 
Hamiltonian which describes a rough surface. 

Above (3e = In (2) the operator cos (2-kK) becomes relevant. Since this is an 
even operator its coupling is nonzero at zero staggered field. This causes the 



model to no longer renormalize to the Gaussian model (as a consequence (4.9) 
is not valid in this region). If (3e > In (2) the surface is thus in a smooth phase 
even if (3s = 0. To complete the phase diagram we must find the behaviour 
of the model for finite values of the staggered field below (3e = \ In (2 — \/2) • 
Before we do that we shall first calculate the singular part of the free energy 
above (3e = | In (2 - \/2) at (3s = 0. 

4.2 Singular part of the free energy of the stag- 
gered F-model 

We sh all a ssume that the staggered F-model can be mapped to a model of the 



form (4T) such that the couplings y„ are analytical as a func tion of (3s in a 



neighborhood of (3s = 0. If the free energy F of the model (4J) is written as 



F = F S +F r (4.14) 
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with F s the singular part of the free energy and F r the regular part of the 



free energy, F s will satisfy the homogeneous part of (4.7) and F r will be a full 
solution of (4.7). Exceptions to this rule may arise when a critical exponent 
associated with the singular behaviour of the free energy is an even integer as 
we shall see later. Ignoring these exceptions for the moment, we see that F s 
satisfies the equation: 



F s (j/i (t) ,y 2 (*)•••)= e2tF s (Vi (0) ,2/2 (0) . . .) 



(4.15) 



In order to see how irrelevant operators modify the singular behaviour, it is 
enough to keep just one irrelevant coupling. The generalization to more irrele- 
vant couplings is trivial. Suppose that for /3s sw the staggered F-model model 
is mapped to a model (4.1) with yx (0) the coupling of the relevant operator 
cos (irh) and y 2 (0) the coupling of an irrelevant even operator. The mapping to 
the model (4.1) can then be written 



yi (0) = Rxf3s + O ((/3s) 3 ) 
tfn(0) = Rn + O ((/3s) 2 ) 



The renormalization group equations (4.2) can be rewritten as 



dyi 
dt 



aiVi 



(4.16) 



(4.17) 



With a\ = 

group equations have been ignored. From (4.15) and (4.17) it then follows that 



-r= > and a n < 0. Higher order terms in the renormalization 



F s ( yi (0) , y 2 (0)) = e~ 2t F s (yx (0) e Qlt , y 2 (0) e a2 *) 
Now choose t such that 

yi(0)e ait = c 



(4.18) 
(4.19) 



with c a constant ^ 0. we can then rewrite (4.18) as 

'yi(oy 



^(1/1(0), 2/ 2 (0)) 



F s [c,y 2 (0) 



Vi 



(4.20) 



Since we expect F s to be analytical as a function of y 2 (0) as long as yx (0) ^ 0, 
we can expand: 



F, 



c,J/ 2 (0)(!) ^]=A + By 2 (0){^) •' +.. 



(4.21) 



Inserting this into ( 4.20 ) and using ( 4.16 ) yields for the leading singularity in 
the free energy (i*\ (/3s)): 

Fi(/fo)~|/fr|£ (4.22) 

while the irrelevant operator contributes a singularity (F 2 (/3s)) of the form 

F 2 (/3s) ~ Fx (/3s) \/3s\^r (4.23) 
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Note that a\ = 2 — -p, and j is given by (4.9). As /3e approaches | In (2 — \/2) 
from above a\ tends to zero, and the singularity in the free energy becomes 
weaker and weaker. What happens at /3e — \ In (2 — \/2) and below is the 
subject of section 4.3 The above argument can easily be generalised to take 



account of the prese nce o f more irrelevant operators and highe r ord er terms 
in the iden tifica tion ( 4.1q ), the renormalization group equations ( 4.17 ) and the 
expansion ( 4.21 ) . By applying a general result |36|] to this case, we find that the 
free energy contains singularities of the form 



F s (Ps) ~ |/3sf 



D + — — > - JS ^ J - 



(4.24) 



where the rii are positive integers. If the exponent becomes an even integer we 
have to multiply the r.h.s. of (4.24) with In \(3s\. We can demonstrate this in the 
case of the leading singularity as follows: According to (4.7) the renormalization 
group equation for the free energy is given by 



dF 

~dt 



2F ■ 



-c(z/i) 



(4.25) 



c is an even analytical function of y\ , because cos (7r/i) is an odd operator while 
the constant operator is even. We are now assuming that — = 2n where n is 



an integer. (4.17) gives 

yi(*)=yi(0)e ait 
If — c contains a term Kyf n , ( 4.25|) can be rewritten as 



dF 

~dt 



2F = K( yi (0)f n e 2t 



From this equation it follows that 

F(y(t)) = K( yi (0)) 2n te 2t +F(y(0)) 



Using ( [4.16| ) and ( [4.26| ) it then follows that 

F s ~f3s 2n \n\[3 S \ 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



It is interesting to see what happens if we let — approach the value 2n. If 
we put — = 2n — e for small e, we can rewrite (4.27) as: 



dF 

— -2F: 

dt 



K( yi (0)f n e^ +ea ^ 



expanding the r.h.s. of this equation in powers of e gives 

^-2F = K( yi (0)) 2n e 2t [l + e ai t + 
dt 

And the singularity in the free energy can be written as 
F s (/3s) ~ I3s 2n [in \0s\ + J In 2 \/3s\ + . . 



(4.30) 



(4.31) 



(4.32) 
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4.3 The case (3e < \ In (2 - V2 

To complete the phase diagram we must obtain the behaviour of the model at 
finite values for the staggered field. This requires us to study the effects of 
higher order terms in the renormalization group equations. To find the most 
important higher order terms we look for terms that are second order in y\. 
These terms are involved in the generation of operators. The most important of 
these terms is the one involved in the generation of the most relevant operator. 
We also look for the lowest order term in the generation of yi arising from 
interactions with operators which have as low an scale dimension as possible. 
To second order in y\ only even operators are generated and the most relevant 
of these is the operator Oi (h) = a 2 (V/i) — -4. It is also this operator which 
through interaction with the operator 0\ contributes to the generation of 0\, 
which is also a second order effect. Since O2 is Gaussian we can calculate this 
effect simply by perturbing the Gaussian interaction j. Denoting the coupling 
of 0\ by y and the coupling of O2 by — 4-, we can write 

m - (2 _ *_ju (4 - 33) 

Although the function A (j) can be calculated using the methods developed 
in the previous chapter, for our purpose we can afford to leave this function 
undetermined. At j = j, which corresponds to (3e < \ In (2 — \/2) , the operator 
0\ is marginal (i.e. right on the boundary between relevant and irrelevant). To 



investigate the phase diagram around this point, we put j = ^ in (4.33) and 
expand in powers of f. To leading order we find 



&■ Ay 2 

iy _ 16 , 



\J, (4.34) 

at 7T j y 



where A = A (fj. Note that these renormalization group equations are similar 
to those for the XY model (see [[[^[Fi|). To be able to construct the phase 
diagram, we must know how to relate j' and y to the model parameters /3e and 



/3s of the staggered F-model in a nonzero staggered field. According to ( 4.16 ) 
we can put 

2/(0) = R((3e)f3s + o((f3s) 3 ) 

I V 2 \ ' 4 - 35 

j'(0) = R'(f3e)+0[((3s) 2 ) 

where we have used the fact that O2 is an even operator. The function R' (f3e) 
can be calculated by using the fact that at zero staggered field the model renor- 



malizes to a Gaussian model with coupling j given by (4.9). We thus find that 



R' (J3 e ) = - arccos ( 1 - -e 2/3e ) - (4.36) 



We now put (3e = \ In (2 — v2j — u in ( 4.16 ) and expand to leading order. We 
find 

2/(0) = Bfis 

f(0) = -(V2-l)u (4 - 3?) 
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where R = R (/3e = ~ In (2 - \/2))- According to ( |4.34| ) it follows that if (i), 
defined as 

1G 



K(t)=y(ty 



ttA 



f ar 



(4.38) 



is a conserved quantity under renormalization. Above j — all flow lines, 
irrespective of the value of K, renormalize to infinity. Below j' = the situation 
is different. Flow lines with negative K end up on the Gaussian line, while flow 
lines with positive K renormalize toward infinity. The flow lines with K = 
thus mark the boundary between the rough phase and the smooth phase below 
j' = 0. Using ( f.37| ) and (4.38), we see that the lines 

j3s = ± - ^__ f-y/2 -l)u (4.39) 



R-JixA 



with u > are the critical lines of the staggered F-model. Points chosen between 
these lines renormalize toward the Gaussian line, points outside this region will 
not. 

We now proceed with a derivation the singular part of the free energy. As the 
critical line is approached from the smooth side, we expect singular behaviour 
of the free energy (note that points on the critical lines itself renormalize to 
the point j' = on the Gaussian line, there is thus no singularity when the 
critical line is approached from the rough side). Since all points on the critical 
lines of the staggered F-modcl flow toward the same point on the Gaussian line, 
critical behaviour is the same all along the critical lines. We can thus content 
ourselves with a calculation of the singular behaviour of the free energy at (3e = 
i In (2 — \/2) as we let (3s approach zero. In this case we are again in the area 
where the identification ( [4.37] ) and the renormalization group equations ( 4.34 ) 
are valid. The initial values are thus j' (0) = and y (0) = Rf3s. According to 
( |4.38 ) we find that K = R 2 {(3s) for the streamlin e tha t passes through this 
point. Eliminating y (t) in favour of j' (t) and using (4.34) gives us the equation 



— = AK + — j 
at ir 

This differential equation is easily integrated: 



t = 



4V ' A-kK 



■. arctan 



\Va^k 



fit) 



Using ( |4.39| ) and the fact that K = R 2 (/3s) z , we can rewrite (|4.41 
tt(1 + v / 2)tan(6») 



t 



16\(3s\ 



" rctan 1+ m n{d) ] ' {t) 



(4.40) 



(4.41) 



(4.42) 



where 9 is the angle at which the critical line intersects the line [3s — 0. Applying 
( 4.18 ) to our case yields the leading singularity in the free energy F s : 



2 (l + V2) tan(fl) 



F s (f3s) ~ e 1810s 



(4.43) 



We can see that the singularity is of infinite order. This is characteristic of the 
Kosterlitz-Thouless transition. Numerical studies using transfer matrix tech- 
niques have yielded similar results on the phase diagram of the staggered F- 
model fl2l. 
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Note that all results have been obtained by using the information present 
in the behaviour of the height-height correlation function of the F-modcl. To 
obtain more results we clearly need more information. In the next chapter we 
shall discuss a simple method that allows one to expand the free energy about 
the line (3e = \ In (2). This expansion can be used to generate more information 
about the mapping of the staggered F-model to Gaussian models. 
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Chapter 5 

Expansions about free 
fermion models 



In this final chapter we will first present Baxters solution of the staggered F- 
model on the free fermion line (i.e. the line (3e = ^ln(2)). Then we proceed 
by expanding the free energy of the staggered F-model about the free fermion 
line. We shall obtain an explicit expression for the free energy to first order. 
By comparing the singular behaviour of this expression to that obtained from 
renormalization group arguments, we are able to verify the known behaviour 
of the Gaussian coupling to first order about (3e = \ In (2) . To simplify the 
computations of the higher order terms we derive a linked cluster method. 

5.1 Definition of free fermion models 

Baxter has solved the staggered F-model at the temperature /3e = \ ln(2) ||. 
Later it was found that this solution could be generalized to other models if 
a certain condition concerning the vertex weights is met. This condition is 
called the free fermion condition because for eight-vertex models satisfying this 
condition the problem leads to a problem of noninteracting fermions in the S- 
matrix formulation. Let u>i be the vertex weight for a vertex of type i (see fig. 



5.1), then the free fermion condition for six-vertex models is: 

W1W2 + W3W4 — w$wq = (5-1) 

The weights Wi may be chosen inhomogeneous. We now proceed by presenting 
a simplified version of Baxter's solution of the staggered F-model. 

5.2 Baxter's solution of the staggered F-model 

Divide the lattice into two sublattices A and B. Choose the vertex energies 



as indicated in fig. 5.1. Consider the ground state in which all A vertices are 
vertices of type 6, and all B-vertices are of type 5. Any state can now be 
represented by drawing lines on the lattice where the arrows point oppositely 
to the ground state configuration. In terms of these lines the six vertices are 
represented by vertices with either no lines, two lines at right angles, or four 
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Figure 5.1: The six vertices and their energies. The upper and lower signs correspond 
to sublattice A respectively B. 



lines. The energies of these vertices are respectively — s, e and s. The next 
step is to replace the original lattice by a decorated lattice by repl acing each 
original vertex by a "city" of four internally connected points (see fig. 5.2). The 
lines on the original lattice are regarded as dimers on the external bonds of the 
decorated lattice. For any configuration on the original lattice, it is possible to 
place dimers on the internal bonds of the decorated lattice, so that the lattice 
becom es co mpletely covered. Now associate to each dimer a weight as indicated 
in fig. 5.2. We now have to choose these weights such that a close-packed 



dimer problem formulated on the decorated lattice is equivalent to our original 
problem. It is a simple matter to see that for this to be the case, we can put 



C = D = E = F = 



0e 



- c -h0s 



iiB(2) 



(5.2) 
(5.3) 

(5.4) 



Note that (5.4) is indeed consistent with the free fermion condition (5.1) 



5.3 The Pfaffian method 

To solve the close-packed dimer problem, we use the Pfaffian method pft EMBg] . 
This method is applicable whenever the lattice is planar, and works by express- 
ing the partition function of the problem as the square root of the determinant 
of an antisymmetric matrix (a Pfaffian) . 

A contribution to Z 2 can be written as the product of two dimer coverings 
C and C . If C connects a point i with a point j, we write 



C(i)=j 



(5.5) 



It is clear that this defines a bijective map on the lattice. C and C" divide the 
lattice into disjoint loops and pairs of neighboring points (bonds) as follows: If 



C(»i) = 


= *2 


C'(k) - 


= h 


C(»3) = 


= u 


C" (u) -- 


= ib 


C(i 5 ) 


= H 


C(* n _i) = 


— 1>n 


C (i n ) ■- 


= i\ 



(5.6) 
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then the points i\ . . . i n form a loop. If n = 2, we don't get a loop but instead 
a single bond. Note that n is always even (even on lattice types on which loops 
containing an odd number of points exist, the loops generated by C and C" 
always contain an even number of points). Since for each loop one has two 
choices to define the actions of C and C within the loop, a given partition of 
the lattice in loops and bonds is consistent with many different configurations 
C and C . Z 2 can thus be calculated by summing over all partitions of the 
lattice in loops and bonds. The contribution a partition makes is given by the 
appropriate product of the weights of dimers, multiplied by a factor 2 L , where L 
is the number of loops in the partition. If we orient each loop, and sum over all 
oriented loops, the factor two for each loop can be ommitted. Now a partition 
of the lattice in oriented loops and bonds defines a permutation of the lattice 
points. For arbitrary points i and j on the lattice, we define Wij as 



Wij — if i and j are not connected, 

W»j — weight of the dimer connecting i and j 



(5.7) 



In terms of the matrix W, we can write: 

z2 = zZ'U W ^U) (5-8) 

*■ i 

where the sum is over all permutations that contain only cycles of even lengths 
(this restriction is denoted by the prime) and the product is over all lattice 
points. Note that the restriction on the summation is only necessary for lattices 



where loops of odd lengths exist. We now want to rewrite the r.h.s. of (5.8) as 
the determinant of a matrix. This is possible if the lattice is planar, and works 
as follows: One tries to factorize the missing sign of the permutation 7r (s (tt)) 



in the sum in (5.8), so that we have 

s M = II s J>t?) ( 5 - 9 ) 

3 

with the Si j depending only on i and j, and Sij = ±1 (we only need to define 
the Si j when i and j are connected). We shall see that a proper choice of the 
Sij allows one to lift the constraint in the summation in ( 5. J). Anticipating 
this result we can write: 

Z 2 = dctR (5.10) 

where 

Ri,j = SijWij (5.11) 



The Sij have to be chosen such that (5.9) is valid for all permutations making 



a nonzero contribution to ( p.8[ ) . The cycles of such a permutation are precisely 
the oriented loops of even length and bonds, and they all have a sign of —1. We 
thus try to define the Sij such that for a closed loop of even length or a bond 
consisting of the points i\ . . . i n we have 

n 

n **.<*+i = - 1 ( 5j2 ) 
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where i n +i — i\. The case n = 2 yields 

S i,3 = ~ s j,i (5.13) 

so that R is antisymmetric. We can now see that permutations containing cycles 
of odd lengths make no net contribution to det R because reversing such a cycle 
changes the sign of the contribution. A permutation that contains a cycle with 
an odd number of points in its interior also makes no net contribution, because, 
the lattice being planar, these points are permuted amongst themselves, so 
that the permutation contains at least one cycle of odd length. We thus have 



to satisfy (5.12) only for loops with an even number of points in its interior. 



This is fortunate, because it is impossible to choose the Sij such that (5.12) is 
satisfied for all loops of even lengths. The Si.j can, however, be chosen to satisfy 
the condition: 

n 

\{si k ,i k+l = {-l) r+l (5.14) 

fe=i 

where r is the number of points inside the loop. It is clear that if the Sij satisfy 
( |5.14| ) for all loops we indeed have Z 2 — det R. We now specialize to the case of 
the staggered F-model. In this case we are dealing with the lattice shown in fig. 



5.2. Choosing the Sij amounts to giving each bond an orientation so that St 



positive if i points to j . The arrows drawn on the bonds in fig. 5.2 represent such 



an orientation. We will now proof that this choice of the orientations satisfies 



the condition (5.14). The proof proceeds by induction, and depends on the 
fact that loops sharing part of their boundaries may be combined to produce 
larger loops. Note that a loop can be broken down into smaller loops if and 



only if the loop has bonds in its interior. On the lattice (5.2), there are two 
types of loops that cannot be broken down. These are the loops formed by four 
internal bonds of a city, and loops connecting four cities formed by four external 



bonds and four internal bonds. For these loops it is easily verified that Q5.14) 



is true. Since any loo p can be broken down into loops of the above type, we 



have to proof that if ( 5.14 ) holds for two arbitrary loops sharing part of their 
boundaries, it also holds for the combined loop. To see this, suppose that there 
are two loops (L\ and L?) with respectively n and ri interior points, with a 
continuous common boundary consisting of q points. The combine d loop (L3) 



will then have r 3 = r\ + r-i + q — 2 interior points. The product in ( 5.14 ) for a 
loop Li will be denoted as s (Li). We then have 

S (L 1 )s(L 2 ) = s(L 3 )(-l) q - 1 (5.15) 

because if we travers L\ and Li in the same direction, we travers all the bonds 
of L3, while the q — 1 bon ds on the common boundary are all travers ed fr om 
both directions. Assuming ( 5.14 ) holds for L\ and L2, it follows from ( |5.15|) : 

a (L 3 ) = (_i)''i+^+«- 1 = (-1P+ 1 (5.I6) 



5.4 Calculation of the free energy 

We have seen that solving the staggered F-model at /3e = 5 In (2) reduces to the 
evaluation of the determinant of the matrix R. To set up a perturbation theory 
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Figure 5.2: The "cities" on the decorated lattice. A and B refer to the two sublattices. 
The meaning of the orientations on the bonds is explained in the text. 



about j3e = \ In (2), we also need the inverse of R. Both the determinant and 
the inverse of R are easily calculated using the following procedure: At each 
vertex i on the decorated lattice, we associate a variables Xi and x\. For x' i a 
given set of variables, we attempt to solve the equation: 



R 



t,]^] 



(5.17) 



We can rewrite this as follows: Introduce coordinates (n, m) on the original 
lattice, so that increasing n (m) corresponds to moving to the right (upward). To 
each city on the decorated lattice we assign the coordinates of the corresponding 
vertex of the original lattice. The variables Xi and x\ are now given by placing 
variables a n>m , 6„, m , c„, m , d„, m respectively a^ m , b' nml c' nm , d' nm for every n 
and m on the four points of the city with coordinates (n, m) as indicated in fig. 
5.3. (5.17) thus becomes 



^<3,n,m 


+ 


W>Vn,m 


- 


*-'G"n— l,m 


— C n,m 


UQ>n,m 


+ 


^^11,171 


+ 


^C-n-\-l,m 




^CVrn + l 


- 


^^"n^m 


+ 


W>U j n,m 


— (1 n , m 


~(-j&n,rn—l 


- 


u>Cn -m 


- 


^^"11,171 


- V 



(5.18) 



We now perform a Fourier transformation on the variables: 



£« 



,-2**(tP+t5P) 



(5.19) 



The Fourier transform of the other variables is defined similarly. In terms of the 



Fourier transformed variables ( 5.18 ) reads: 
ua v . a + ub v „ - 



Ca% Piq 
-Ca~ q a I} 



ub p 
uc p 
uc v 



Cp- p d p , q 


— C p. q 


Cf3 p c p , q 


— u, pq 


UCLp^q 


— ft p q 


IKXp.q 


— b' 

— U p,q 



(5.20) 
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Figure 5.3: The variables a n ,mi b n , m , c n ,m and d n ,m, associated with a city on the 
decorated lattice. 



Here a = e M and [3 = e « . The determinant A Pi? of ( 5.2C ) is given by: 



A M = 2 cosh (2 j3s) + 2 cos (^ j cos ( ^ 



(5.21) 



And the reduced free energy per vertex (i.e. the free energy times — 0, denoted 
as F) for an infinite by infinite lattice follows: 



F = limjv 



M— »oo 



7X7 In det R = lim 



27V M 



N,M^<x> 



JU Z^p,g m ^V< 



2NM 



giy /„ /„ In [2 cosh (20s) + 2 cos (0i) cos (0 2 )} d0id0 2 



(5.22) 



5.5 Singular part of the free energy 

In this section we calculate the singularity in the free energy of the staggered 
F-model at (3s — on the free fcrmion line. We shall use the following method: 
We expand the logarithm in (|5.22), thereby obtaining a series of the form 



F(Ps) = J2 



n 2 n 3 



^ cosh 2 " (20s) 
Summing the series inside the brackets term by term we obtain 



(5.23) 



F(/3s)=]TF fc (/3 s ) 



k=2 



with 



(5.24) 



F k (0s) = £ 



\ n k cosh^ n (20a) 



(5.25) 
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As can be seen by differentiating Fk repeatedly, the singularity in Fk becomes 
weaker as k increases. 

Expanding the logarithm in (|5.22[) yields 



F(f3s) 



Sw 



1 r* {"* rift rift Y^°° C0S"(6> 1 )C0S"(6' 2 ) 

"" J-7T 1-tt uVlUfy 2^„=1 ncosh"(2/3s) 



1 ^—vOO 

2 



Eoo 1 
n=l 2rTcosIr ZT 



(2/3s) 



(2»)! 



4™n! 2 



Using 



= n n e ™v / 27mexp V^ 



£. 



21: 



\k=l 



2k{2k-l)n 2k - 1 



(5.26) 



(5.27) 



where the B r are the Bernoulli numbers, we find 



F 



oo 



4tt j^j n 2 cosh ^n ( 2 /3a) 



1 



1 



1 



4n 32n 2 128n 3 2048n 4 



(5.28) 
We can find the singular part of the function V°° , — - — ,\ n ,„ a ; as follows: 

° -^ ^-^n— 1 nP cosh^ 71 (2/3s) 

Put t = In (cosh (2/3s)) . We then have to find the singular part of the function 
Up (t) with 



M*)=£ 



??.'-' 



(5.29) 



as i -> for p > 2. From (|5.29| ) it follows that 



rff 



t/„ 



(5.30) 



We denote the singular part of U v by U p . It then follows from ( 5.3C ) that 

dUp+i 



-*r = - u * 



For p = 1 the sum in (5.29) is easily evaluated: 

Ui{t) = -ln(l-e -4 ) 
And we see that U\ is given by 

U x = - In (t) 
From ( |5.33 ) and ( 5.31 ) it then follows that 



Up(t) = (-lT 7 ^—^\nll) 



(P-1)1 



Inserting this in ( [5.28] ) gives 



r 1 



1 / £ 2 i 3 

F s (/3s = it + — + 

yy ' 4?r V 8 192 3072 



ln(t) 



(5.31) 
(5.32) 
(5.33) 

(5.34) 
(5.35) 



Where F s (/ 3s) is the singular part of the free energy and t = 2 In (cosh (2/3s)). 
Expanding (5.35) in powers of (3s gives 



F s (/3s) 



<« 2 -> 4 +li <«*-!>> 



\n\f3s\ 

(5.36) 
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5.6 Perturbation theory 

We now proceed with the derivation of a perturbation theory about the free 
fermion line of a 6-vertex model. The Hamiltonian of a general 6-vertex model 
can be defined as follows. One assigns an energy e(p,i) to a vertex in state p 



(see fig. p.l[ ) and position i. The configuration of the lattice can be specified 
by a function c which maps a position of a vertex to a number, f • • • 6, which 
is to be interpreted as the state of the vertex at that position. The reduced 
Hamiltonian (H) is defined to be the functional that assigns to each state c its 
energy times —f3. We can thus write 

fl r (c) = -/3^e(c (»),*) (5.37) 

i 

For H a Hamiltonian of a general 6-vertex model and H a Hamiltonian of 
a free fermion model a perturbation V can be defined so that we have 

H = H + V (5.38) 

The partition function Z can be written as: 

Z = J2 eH ° ic)+V(c)=Z o(e V ) (5-39) 

c 

Here Zq is the partition function of the free fermion model. The reduced free 
energy can be expressed as: 

F = F a + In (e v ) = F a + (V) + i ((V - (V)) 2 ) + ... (5.40) 

Here Fq is the reduced free energy of the free fermion model. Now write V — 
J~). Vj with Vi (c (i)) a perturbation of the vertex energy times — /3 at position i. 
( |5.40| ) can be rewritten as: 

F = F + J2 W> + \zZ WW - W ^>] + ■ ' • ( 5 - 41 ) 

i vj 

To compute a free fermion average (V^V^ ...V% n ), we can proceed as follows: 
Introduce a constraint in the free fermion model by requiring the vertices at the 
positions i\ . . . i n to be in the states x\ . . . x n . The partition function of this 
model is denoted by Zt t ...i n (x\ . . . x n ). We can then write 



(V h V i2 ...V in )= J2 



Z ix ... in (xi . . . x n ) V (aji) ...V (x n ) 



x 1 ...x n ^ (5.42) 

It now remains to calculate Zi 1 ...i n (x\ . . .x n ). It is convenient to reformulate 
this problem as follows: Denote the state of an arrow located at the bond j by 
Sj . Put Sj = 1 if the arrow points oppositely to the ground state configuration 
and Sj — otherwise. Define a constrained free fermion model by requiring the 
arrow at the bond j r to be in state Sj r for 1 < r < m. We then want to evaluate 
the partition function of this model, which we denote as Z (sj 1 . . . Sj m ). The 
idea is to perturb the weights of the dimers on the bonds j r infinitesimally. We 
assign a weight C (1 + e r ) to the dimer on the bond j r . The partition function 
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of the (unconstrained) free fermion model [Z (ei . . . e m )) can be written in terms 
of the constrained partition functions as: 

Z(ei... e m ) = J2 {s } z ( s h ■ ■ ■ Sj m ) rifcli (! + s j^k) 

= z o + J2k Z ( s j k = !) £ fe + Efc<z z ( s ik = h s n = !) £ fc£/ + ■ ■ ■ (- 5 43 ) 



Z (ei . . . e m ) can be calculated using ( 5.10 ), by making the necessary changes to 
R. We can write: 

n i 

R = R + J2^kR(k) (5-44) 

fc=i 

Here i?o is the original unperturbed matrix, Rqa is dehned as follows: 

R(k),ij = c 



if i and j are connected by ju and i points to j, 

R(k),ij = -C 

if z and j are connected by jk and j points to i, 

R(k),ij = 
if i and j are not connected by jk . 



Note that the Rtk) have only two nonzero matrix elements. Inserting (5.44) in 
(|5.10) and expanding gives: 



z = Vdetl != % /d^n^ Trln [ 1 +E^ fl o- 1 ^>] 

= Vdetl^ [1 + | E* e * Tr (V#(fc)) 

+1 £*,, c*Q [5 Tr (iio" 1 ^) Tr (i^,)) - Tr (i^ 1 ^^ 1 ^))] + .(.5,fe) 



Using (5.45J) and (|5.43) we can directly read off the constrained partition func- 



tions which have all their arguments set to +1 (i.e. all the constrained arrows 
point oppositely to the ground state configuration). To calculate a constrained 
partition function with some of its arguments set to 0, we simply have to apply 
the principle of inclusion and exclusion (i.e. a Mobius inversion on the power 
set of arguments). For example consider the evaluation of Z (si, S2, S3, S4, S5), 
with Si = 82 = 1 and S3 = S4 = S5 = 0. Put £3 = £4 = £5 = 1. According to the 
principle of incusion and exclusion, we can write: 

Z(si,S2,S3,S 4 ,S 5 ) = Z(si,s 2 ) 

- [Z(si,s 2 ,t 3 ) + Z (si,s 2l t 4 ) + Z(si,s 2 ,t 5 )] 

+Z (si, s 2 , t 3 , £ 4 ) + Z (si,s 2 ,t 3 ,t 5 ) + Z (si, s 2 , £4, £5) 

-Z(s 1 ,s 2 ,t 3 ,t 4 ,t 5 ) (5.46) 

5.7 First order computation for the staggered 
F-model 

For the staggered F-model the expansion can be simplified. The vertex in the 
ground state at a particular point will be referred to as an a-vertex. A b-vertex 
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is obtained by reversing the arrows of an a-vertex. An a-vertex (b-vertex) is 
thus of type 5 or 6 and has an energy of — s (s). The constrained partition 
function corresponding to the model with one vertex constrained to be an a- 
vertex (b-vertex) is denoted as Z a (Zb). Note that under the transformation 
S — > — s the role of vertices a and b are interchanged. We thus have 



Z a (/3s) = Z b (-/3s) 



(5.47) 



If we put fit = ^\n(2) + U we have, according to (5.41) and fl5.42 ), to first order 
in U: 



F = F Q 



z ~ z :~ Zb u+o(u*) 



(5.48) 



Here F is the reduced free energy per vertex of the staggered F-model. To 
calculate Zb we only have to constrain two opposing arrows of one vertex to 
point oppositely to an a-vertex. If we choose to constrain two horizontal arrows, 
we need the i?~ for i and j corresponding to an a or b variable on the decorated 



lattice. We define Green's functions G 



(*>,x' { 



j) (p, q) with a;W € {a, b, c, d} and 



>U) 



6 {a', 6', c', d'}, so that the solution of (|5.20) can be written as 



,W 



XXwy^M^' 



C?) 



(5.49) 



Solving equation (5.2C) yields 



G a y (p, q) = 




- 


1 e ^ l3a (8 P 3~ p ) 

A(p,q) 2 yl J P > 


G a ,b> (p, q) = 


i 

A(p,<?) 


4^ (p + /3~ p ) + e~^'ofl 


Gb, a > (p, q) = 


1 
A(p,g) 


^([3P + f3-P) + e-^ s a-i 


Gb,b' (p, q) = 






1 a 

1 e5" lav o-v\ 
\(v.a) 2 \P P > 



(5.50) 



An inverse Fourier transformation yields in the limit N,M — ► oo: 



G a ,a> (n,m) 
G a ,b>{n,m) 

Gb, a ' (n,m) 



hr Jo* Jo* fflifflv 



,( e .(»+l)8 



.(n-l)fl 2 1 



&S?j? r d0 1 d0, 



A(0i,0 2 ) 



^! e ™»l( e .(» + l)«2 +e . 



(n-l)fl 2 '\ 



-££* .(m+l)« 



i^ft** 



A(0 1: 2 ) 



4£f. e .ms 1 ( e .(»+i)« 2+e 

A(0i,0 2 ) 



.(n-l)S 2 N 



(5.51) 



G b ,b'(n,m) = %j3- /„" | * d(M6>2 

Here A (#i, 2 ) is given by 



A(0 1: 2 ) 



A( 



A (6»i , 6» 2 ) = 2 cosh (2/3s) + 2 cos (6>i ) cos (6 2 ) 



(5.52) 
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The necessary components of the matrix R x can be expressed in terms of the 
Green's function G as: 



R i = G x x ' (n — n! ,m — m!) 

Put W h = J|. It follows from (|1|) and (§1||) that 

W b = ~ Tr 2 (i^ 1 ^) - ~ TV (^R^^R. 



(')) 



(5.53) 



(5.54) 



where fc and / refer to the bonds that connect the points a n>m and 6ra,m+i 
respectively a n , m +i and b n , m +2 on the decorated lattice. Using (5.7) the first 
term in (5.54) (lUbi) can be written as 



W b i = e —^ [G bta (0, 1) - G a . b (0, -l)] 2 



(5.55) 



Inserting the appropriate expressions given in ( 5.51 ) in this equation yields 

I / -" / i ■ "■ A- en-, (f/, i ens I H-, I 

w bl 



64tt 4 



i, .;,, cosh (2/3s) + cos (0i) cos (#2) 



(5.56) 



The second contribution to Z b (W&2) can be expressed in terms of the Green's 
function G as: 



W b2 



■s-f- [G ba (0, 2) G ba (0, 0) + G ab (0, 0) G ab (0, -2) 



G 



aa (0, 1) G b6 (0, -1) - G aa (0, -1) G bb (0, 1)] (5 5?) 



Using ( 5.51 ) it is not difficult to see that each of the four terms in the brackets 
of this equation vanishes. We thus have 



W b = W bl 



1 



64tt 4 
From ( |5.47| ) it follows that 



27 '■-* e" 2 ' 3s + cos(6'i)cos(6' 2 ) 



,, ./,, cosh(2/3s) + cos (0i ) cos (02 ) 



W a (13s) = W b (-0a) 



(5.58) 



(5.59) 



The expression obtained by substituting ( 5.58J ) and (5.59) in (5.48) can be sim- 
plified by exploiting the fact that we only need the even part of the function 
Wb ((3 s). We will denote even (odd) parts of functions with a subscript + (— ). 
( |5.48|) can be rewritten as 



F = F Q + (2W b+ -1)U + 



( [5.58 ) can be rewritten as 



W b = I 2 



with 



e- 2/3s + cos(6>i)cos(6>2) 



rlf) HO 
,, ./,, cosh (2 f3s) + cos (81) cos (#2) 



(5.60) 
(5.61) 

(5.62) 
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From (|5.61) it follows that 



W b+ =11 +11 



2 , jl 



From ( |5.62 ) it easily follows that 



(5.63) 



(5.64) 



and 



(|5.60) can thus be written as 



F = F 



I_ 



2d/3s 



dF n 
dps 



- 1 



U 



(5.65) 



(5.66) 



5.8 Singular behavior in the vicinity of the free 
fermion line 



Equation (5.66) allows us to check the dependence of the Gaussian coupling j 
on f3e as given by (4.9) for /3e ~ ^ln(2). We have seen that renormalization 
group arguments lead to the following leading singularity in the reduced free 
energy: 



with 



If we put 



F s ~(f3s) 2 -*J 



j = - arccos [ 1 — -e 2/3e 
J 2 V 2 



(5.67) 



(5.68) 



fh = - In (2) + U 



(5.69) 



in these equations, we find that 



A(U)(psY 
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— (U + O (U 2 )) In (ps) + — (U 2 + (U 3 )) In 2 \0s\ + 



(5.70) 



If we compare this with (5.36), we find that the amplitude A (U) is given by 



A ^ = W 



(5.71) 



It then follows that the amplitude of the term {[3 s) 2 In 2 \/3s\ is A (U + O ( U 2 )) - 
It is now a simple matter to verify this using ( 5.66| ) and ( 5.36 ) . From ( |5.36 ) 
and ( 5.66 ) it follows that the order U contribution to the singular part of the 
reduced free energy F\ ((3s) can be written as 



F l<8 (/3s) = [Bi {ps) In \Ps\ + B 2 (ps) In 2 \ps\] U 



(5.72) 
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with B\ and B 2 regular functions of (3s. Inserting (5.36J) in (^>.66 ) gives 



B 2 (Ps) = — 



m 2 ~l(f3s) 4 + 7 ^(/3s) 6 + 



(5.73) 



We have thus verified ( 5.6q) to first order in U. 



5.9 Linked cluster expansion 

In this section we will derive a diagrammatic (linked cluster) method to compute 
terms in the perturbative expansion of the free energy. Instead of perturbing 
the energies of vertices 1...4, we will perturb the energy of an a- vertex. We 
will consider the model with the following vertex energies: vertices 1...4 are 
assigned an energy of eo, with (3eo = ■%- In (2) , a- vertices are assigned an energy 
of — so + v, and b- vertices are assigned an energy of sq. The reduced free energy 
of this model is denoted as F (—(3v). It is not difficult to see that the reduced 
free energy of the staggered F-model can be obtained from F: 



F 



pe = (3e + ^U,f3s = f3s + ^U) =F{U)-^U 



We now define a Hamiltonian (H ({J})) as 



(5.74) 



(5.75) 



Ho is a free fermion model and the Vi are defined as follows. If c is an arrow 
configuration then V, (c) is the increase in energy times — (3 of the vertex at 
position i (cj). In this case we thus have Vi (c) = U if c, is an a-vertex and 
Vi (c) = otherwise. If we put all the J, = 1, H ({J}) becomes the Hamiltonian 
of the model defined above. The partition function for H ({J}) is denoted as 
Z ({ J}). In the following wc shall use the notation J = p to indicate that Ji = p 
for all i. An expansion about J = yields 



Z(J=1) 



oo 1 



d k Z({J}) 



, k\ .^-f dJ lt ■■■dJ ik 

k=0 ii---ik 



(5.76) 



j=o 



From the definition of Z ({J}) it follows that the free fermion correlations 
g (ii • • 'ik)i defined as 



can be expressed as 



'(«1,«2 



g{h,i2- ■■ik)= (Yi t V i2 ■ ■■V ln ) 
d k Z({J}) 



1 



ik) 



z(j = o) dj n ■ ■ ■ dj z , 



(5.77) 



(5.78) 



j=o 



The terms in the expansion ( 5.76 ) are thus readily expressed in terms of the 
free fermion correlations. We have 



OO -. 



(5.79) 



k=0 
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To compute the free fermion correlations we define for i\- ■ -ik k points on the 
lattice, a constrained free fermion model by demanding that the vertices at 
positions i\- ■ ■ ik be a-vertices. The partition function of this model is denoted 
as Zi t ...i k . In terms of Zi 1 ...i k , g(ii ■ ■ ■ ik) can be written as 

»<•■■-■•*> = zttSV 1 ' < 5 - 80 > 

From the correspondence of arrow configurations to dimer configurations, it 
follows that Zi 1 ...i k can be written as the sum of all dimer configurations such 
that at the cities ii • ■ ■ ik dimers are placed on two internal bonds. For each of 
these cities there are two choices for their internal dimer configurations. The 
sum of all dimer configurations consistent with one particular choice at each 
city is denoted as Z' ii ... ik . Since Z' ii ... ik does not depend on the particular choice 
made, it is clear that we have 

Z n ... lk = rZ[ % ... lk (5.81) 

where r is the number of different indices in i\ ■ ■ ■ ik- For definiteness we assume 
that at each city i p , for p — 1 ■ ■ • k dimers are placed on respectively the bond 
connecting the points labeled by the variables a% and a , and di and bi (see 



fig. 5.3). These bonds will be referred to as modified bonds, and the position 
of these two bonds will be denoted as, respectively, i Pl and i P2 . Dimers on 
modified bonds will be referred to as modified dimers, and the modified lattice 
is defined to be the decorated lattice, but with only the modified bonds. To 
compute Z% x ...i k , we define an inhomogeneous free fermion model by changing 
the weights of the modified bonds. At all positions j we associate the variables 
Cjj and ej 2 to, respectively, the modified bonds at j\ and j2- We now choose 
the weight of a dimer on a modified bond at a position A: to be a function of Ck, 
so that we have 

Z ix ... ik =cocff. of e tli e il2 ■ ■ ■ e iki e ik2 m Z ({e}) (5.82) 

Here Z ({e}) is the partition function of the inhomogeneous free fermion model. 
We denote the weight of the dimer on a modified bond at a position k as w(ek)- 
The choice w (efe) — u (l + v2efc) (u is the original weight of the dimer) seems 



to work, but there is a problem if some of the indices in (5.82) coincide. To see 
this, consider the r.h.s. of (5.82|), when all the indices in this equation are chosen 



differently. It is clear that a dimer configuration contributing to Z ({e}) has a 
weight proportional to e^ e^ • * • e» fc £i k only if it consists of all the modified 



dimers at the positions i\ 1 ,i\ 2 ...ik 11 ik 2 - The r.h.s. of (5.82) is thus equal to 



2 k Z' i i , and this equals according to (5.81) the r.h.s. of (5.82). It is also 



clear that with the present choice of the dimer weights Z ({e}) is a multilinear 



function of the e's, so (5.82) is certainly not true when not all of the indices are 
different. Although this doesn't seem to be a big problem, because there is no 
point in constraining a certain vertex more than once, we have to sum over all 
indices to calculate Z ( J = 1) and it is convenient to do so unrestrictively. The 
solution to this problem is obvious. We define w (e&) as 

w(e k )=u(l + ^-\ (5.83) 
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Since Z ({e}) is a free fermion model, we can use (5.10): 



Z({e}) = ^/detR({e}) (5.84) 

with 

R({e})=R(e = 0) + J2R(k) (5.85) 

k 

where the sum is over all modified bonds, e = means e, = for all i and Rqa 
is defined as follows. If the bond k connects the points i to j we put 

R {m = ±V2-^±- (5.86) 

t £k 

The plus sign is chosen if the orientation on the bond is such that i points to j, 
else the minus sign is chosen. All other matrix elements of R^) ar e zero. Q5.84 ) 
can be rewritten as 

Z({e}) = Z e5 Tr Mi+GB) (5 . 87 ) 

where G = R' 1 (e = 0), Z = Z ({e = 0}) and B = £ fc %)• Note that B v is 
only nonzero if i and j are connected by a modified bond, in which case 

Bij = R(k)i 3 (5.88) 

if the modified bond k connects i to j. Expanding the logarithm in this equation 
allows us to rewrite the exponent of the above equation as 



- Tr In (1 + GB) = -J2 ^— £ II G ir , jr . +1 B jr+uir+1 

n=l iiji-in,j n r=l (5.89) 



Here i n +i = i\ and j n +\ = ji- Note that for a contribution to the summation on 
the r.h.s. of this equation to be nonzero, i r +i must be chosen so that j r +i and 
i r +i are connected by a modified bond. The summation can then be interpreted 
as a summation of amplitudes of all closed paths of n steps, where a single 
step consists of going from a position i r to an arbitrary position j r +i, and 
then from j r +i to that point (i r +i) which is connected to j r +i by a modified 
bond. If r = n, j r +i is, of course, not arbitrary because of the conditions 
?'„+! = i\ and j n +i = ji- Hence it follows that we are dealing with closed paths. 
To each step we associate a weight which in the example given above equals 
Gi r j r +iBj r+lt i r+1 . The amplitude of a path is the product of the weights of 
the steps forming the path. Each closed path can be defined as an oriented 
loop of which one point plays the role of starting and end point. It is clear 
that the amplitude of the path does not depend on the choice of the starting 
point. It also does not depend on the orientation, because both G and B are 
antisymmetric. Instead of summing over closed paths, we may thus sum over 
loops. To do this, we need to know the number of closed paths that correspond 
to one loop. It is clear that this is given by two times the number of ways one 
can attach a starting point to the loop. If the loop consists of n steps, and 
if the winding number of the loop is p (i.e. traversing the loop once, means 
visiting the same points in the same order p times) then there are — ways of 
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attaching an starting point to the loop. Suppose all loops are enumerated in 
some arbitrary order. We define the amplitude of the r loop, denoted as L r , 
as the amplitude of one of the closed paths to which it corresponds divided by 
the winding number of the loop and multiplied by a sign (±1). This sign is the 
product of the following factors. For each modified bond there is a factor — 1, 



and there is an overall factor of —1. We can now rewrite (5.89) as 



~Trln(l + GB)=Y^L r (5.90) 

2 r=l 



Inserting this in ( |5.87|) gives 



Z({£}) ^f[e^ ( 5 - 91 ) 



r—l 



It thus follows that if'' is a sum of amplitudes of diagrams, where a diagram 
is a set of loops and the amplitude of a diagram is the product of amplitudes 
of the loops of which it consists multiplied by a factor Yl^Li FT ^ 1°°P r occurs 
k r times in the diagram. Note that with our original definition of the weight 
function w (e) = u (l + v2e) , Z ({e}) would be a multilinear function of the e's. 
This implies that products of loops sharing bonds will not make a net contri- 
bution (this does not, of course, depend on the choice of the weight function). 
To see an example of this consider two different loops with amplitudes L\ and 
L2 sharing a bond. In that case there is a contribution L1L2 to Z ({e}). There 
is also a single loop, obtained by merging the two loops. This loop has an am- 
plitude of —L1L2, and we see that the two contributions have indeed cancelled. 
Similarly, contributions from loops with winding numbers greater than one can- 
cel against products of loops with smaller winding numbers. Although, as we 
shall see later, such cancellations will also occur in the expansion of Z (J = 1), 
ln(Z(J = 1)) will contain only connected diagrams (the precise definition of 
"connected" will be given later ) . Diagrams in which a bond appears more than 
once can make a net contribution to ln(Z(J = 1)) because the disconnected 
diagrams against which such a contribution would cancel don't appear in the 
expansion of In (Z (J = 1)). To take into account the cancellation of loops with 
a winding number greater than one against products of loops with lower winding 



numbers, we can rewrite (5.91) as 



Z{{e}) -f{[l + L r } (5.92) 



Zq 

r—l 



Here the product is only over loops with a winding number of one. To derive 
this equation directly from (5.90|) one proceeds as follows. The amplitude of a 



loop is calculated as above, but without the overall minus sign. We may then 
write the amplitude of a loop with winding number p as 

A r , p = — (5.93) 

Here A T:P is the amplitude of the loop with winding number p an d A r is the 



"same" loop with winding number one. We can then rewrite (5.90) as 



DC 



-Trln(l + GB) = -J2J2f=z2 ln ( 1 - A r) ( 5 ' 94 ) 

r=l p— 1 r—l 
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Here the sum over r is only over loops with a winding number of one. Inserting 
this in (5.87) then gives ( 5.92Q (In that equation L r is defined as —A r ). With 
our definition of the weight function ( |5.83| ) a diagram contributing to Z ({e}) 
is a multinomial in the e's. To have that each diagram is a monomial, we split 
each modified bond on the decorated lattice into an infinite number of bonds, 
each of them connecting the same two points as the original bond. These bonds 
are all assigned the same orientation as the original bond. The weights of these 
bo nds a re chosen by assigning to each bond a (d ifferent) term in the expansion 
of ( 5.8S ) in powers of e. It now follows that ( 5.91 ) still holds on this new lattice. 
We can see this as follows. Let all loops on the new lattice be enumerated. The 
amplitude of the k l loop will be denoted as Lk- Since the sum of the weights 
of the new bonds replacing one modified bond equals the weight of the modified 
bond, we can write 



zJ Lr - Yl ~ Lk 



r=l 



k=\ 



( |5.91 ) can thus be rewritten as 



(5.95) 



Z({e}) 

Z 



n 



.Lr 






n 

fc=l 



(5.96) 



Note that the cancellation of diagrams discussed above implies that a diagram 
on the new lattice which uses more than bond connecting the same points will 
make no net contribution. In particular loops with winding numbers greater 
than one will cancel against products of loops with smaller winding numbers. 
To emphasize this we can rewrite ( [5.92] ) as follows. We say that a loop on the 
new lattice corresponds to a loop on the modified lattice if the same points can 
be traversed in the same order. The amplitude of a loop on the new lattice 
corresponding to a loop r on the modified lattice (both with a winding number 
of one) shall be written as L r ^, where the index k enumerates all such loops. 
We can then write 



Inserting this in ( J5.92j ) then gives 



fe=i 



(5.97) 



Z({e}) 

Z 



n 



fc=i 



r.k 



(5.98) 



We shall use the following notations for the bonds in a city. For a city at 
position i, there is a bond with a weight proportional to e^ and q 5 , if k > 1. 
These two bonds will be referred to as bonds of order k. A bond that connects 
the points Oj and Ci will be referred to as a bond of type 1, a bond that connects 
the points bi and di will be referred to as a bond of type 2. A diagram will 
be called balanced if at every city the sum of the orders of bonds of type 1 in 
the diagram equals the sum of the orders of bonds of type 2. From (5.82) it 
follows that only balanced diagrams contribute to Z^...i n . It is also clear that a 
particular balanced diagram contributes to Zi x ...i n only if at each city the sum 
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of the orders of bonds of each type equals the number of times the city occurs as 
an index in Zi x ...i n . For balanced diagrams the sum of the orders of bonds of a 
certain type at a city will be referred to as the order of the diagram at that city. 
If a city is not part of a diagram, then the order of that diagram at that city 



is zero. From (5.80) it now follows that g (i\ ■ ■ ■ i n ) can be written as the sum 
of amplitudes of balanced diagrams, if the weights associated with bonds are 
redefined as follows. We define the weights associated with the modified bonds, 

by choosing e = VU : so that a bond of order k has a weight of ±\/2 I VU) . 

9 (ii • • • in) can thus be written as a restricted sum of amplitudes of balanced 
diagrams, where the restriction implies that a diagram contributes if at each 
city the order of the diagram equals the number of times the city occurs as an 



argument of g. (5.79) expresses Z (J = 1) as 



Z(J = l) = Z(J = 0)J2y E 9(k---ik) (5.99) 

k=0 ' ii---ik 

By substituting the diagrammatic expansion of g in this equation, we see that 
Z /jZ l is a sum of balanced diagrams. We denote the order of a balanced 
diagram at a city j as rrij. We now focus on a balanced diagram with yV - m,j = k. 



This diagram will appear yr^ — - times in ( 5.79] ) . The contribution this diagram 



makes to Z ('j— l is thus given by its amplitude multiplied by the factor *-r 1 

The fact that this extra factor appears is rather inconvenient. This problem can 
be overcome by again changing the weights associated with the bonds. Since 
diagrams which use more than one bond of a certain type at a city cancel, we 
can absorb a factor — t= in the weight of a bond of order r. The weight of a bond 

Vr! 

of order r now becomes ±V2 A ■ With this definition of the weights of the 

Vr! 

bonds z(j=tl) 1S thus simply the sum of amplitudes of all balanced diagrams. 
Note that the factor -4= is rather arbitrary, we could have chosen ?4p for the 
weight of a bond of order r of one type and j-^y; for the bond of order r of the 
other type, provided that p + q = 1 . This arbitrariness shall lead to a drastic 
simplification for the expansion of the free energy. 

We now proceed to write z \jZ l in a way analogous to the expression ( |5.96 ) 



for z , so that it becomes a simple matter to find the diagrammatic expan- 
sion of In (Z ( J = 1)). We can do this as follows. We define a diagram to be 
connected (irreducible), if it is balanced and cannot be written as a product 
of smaller balanced diagrams (by a smaller diagram we mean a diagram that 
contains only a part of the loops of the original diagram) . Suppose all connected 
diagrams are enumerated in some arbitrary order. The amplitude of the r 
connected diagram is denoted as C r . We can now express z \' tZJ. in terms of 
the C r as 

^-^ = f[[l + C r ] (5.100) 

r— 1 

In this equation the product is over connected diagrams consisting of loops with 
winding numbers of one. The reduced free energy of the staggered F-model can 
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thus be written as 

oo 

F = Fo + ^]ln(l + CV) (5.101) 

r=l 

As noted earlier, there is some freedom in the choice of the weights of the 
bonds. At a certain city, only the product of the weights of the two bonds of 
order r is required to contain the factor A . For the expansion of the free energy 
this means that a diagram that uses a bond of order r of one type also has 
to contain the bond of order r of the other type, otherwise the diagram will 
cancel against other diagrams. We can thus go back to the modified lattice by 
summing over the orders of the bonds. The weight of a modified bond then 
becomes ±-i/2 (exp (u) — 1). The sign is given according to the orientation of 
the bond. The extra minus sign can now be omitted, because any diagram has 
an even number of bonds. 
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Appendix A 

Notations and conventions 



A.l Vectors 

Vectors are denoted like ordinary variables (i.e. without an arrow). The modulus 
of a vector x is denoted as |x|. By x n with n an even integer we mean \x\ . 

A. 2 Summation convention 

Repeated indices are summed over, unless stated otherwise. 

A. 3 Multi-indices 

A multi- index containing n indices is denoted as (n). We use a summation 
convention in the following way: 

^(n)-B(ro) = -Ai 1 ...» n Bi I ...i n (A.l) 

A possible sum over n is always explicitly denoted. Different variables between 
brackets always contain different indices, e.g. 

A( n )B( m ) = A il ... in Bj 1 ...j m (A. 2) 



Identical variables between brackets contain identical indices (as in (A.l)). A 
tensor can have more than one multi-index: 

A(n),(m) = Ail— in,h—jm (A-3) 

By 

^2 A"i).(«a)-(«r) ( A - 4 ) 

contractions 
we mean a sum over all possible contractions of the indices with each other, e.g. 

Yl A («i),(«2) ( A - 5 ) 

contractions 

with m — ri2 — 2 stands for 

A-ukk + Aikij- + Aikki (A.o) 
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A. 4 Derivatives 

For derivatives we sometimes use the notation d(t{ >,• By this we mean a the 

k derivative with respect to the components (k) of b evaluated at the point 
a. The arguments a and b are optional. 



A. 5 Distributions 

For a detailed account of distributions we refer to ||31^J3!q| . Once a well defined 
space of functions S is defined, a space of continuous linear functionals S' on S 
can be defined. We shall define the set S as the set of all functions h(x) that 
can be written as in ( |3.8[ ). Apart from the usual L 2 -norm we define a Frechet 
norm as follows: 

\h\= sup \d (k) h(x)\ (A.7) 

(k),x,k<N 

Here N is some arbitrary number. The set of all continuous linear functionals 
(distributions) with respect to this norm is denoted by S' N . The action of a 
distribution T £ S' N on h 6 S is denoted as Th or as T(h). It is not difficult to 
see that the functional T defined as 

Th=d in) h\ z (A.8) 

for some multi-index (n) is an element of S' N when n < N . The functional I 
defined as 

also defines a distribution as long as I(x) is bounded. A weak-topology on S' N 
is defined by saying that T n -»• T if T n h -> T/i V/i <G 5. 
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Appendix B 

Gaussian correlations 



B.l Correlation function for h^ 



In this appendix we derive the expression ( 3.27| ) for the correlation function for 
the field h^ 2 \ By definition we have: 



G (x) = (/i (2) (x) h^ (o) 



In terms of h this becomes 



V 



(B.l) 



(B.2) 



fci,fc 2 



where the sum is over elements of the set S^ (see section |3.3| for the definition 
of this set). Since in the Gaussian model h(k\) and h (fe) are uncorrelated 
unless fei = ±Ai2, we have 



G ^ = V ^ (\h(k)\ 2 + (h(k)) 



ik-x 



kesw 
Now define real valued variables a (k) and b (k) as follows: 

h{k) =a(k)+ib(k) 
Since h (x) is real we have 

h(-k) = a{k)-ib(k) 



(B.3) 



(B.4) 



(B.5) 



From ( p.4[) and (|B.5D it follows that 

dh{k)dh{-k) = 2da{k)db{k) (B.6) 

The partition function of the Gaussian model can be written as 

-j(k)k 2 (a(k) 2 +b(k) 2 ) _. TT '^ N ( k ) 



Z= H / da(fc)d&(fc)7V(A:)e^W fe2 H fe ) 2+fc « 2 )= [] 



U J H&7) 
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The prime indicates the condition k y > and we have replaced the Gaussian 
coupling j by a fc-dependent coupling j (k). N (k) is an uninteresting function 
coming from the definition of the measure Dh (see ( |3.14| )), and the Jacobian of 
the transformation to the variables a(k) and b(k). It is clear from (B.7) that 



(h(k))' 



vanishes. 



h(k) 



is computed as follows: 



h(k) 



1 dln(Z) 



k 2 dj(k) 



Using this one can write (B.3) as: 



<*(*)=4 E 



Kk)=3 



„ik-x 



jk 2 



V ^ jk 2 
keS( 2 ~> 



(B. 



(B.9) 



(2) 



We can formally rewrite (B.9) using the characteristic function P c of the set 
S^ as: 



G(x) 



d 2 k P^ (k)_ elk . x 



(B.10) 



(2tt) 2 jk 2 

B.2 The height-height correlation function 

In this section we will evaluate the height-height correlation function (R(x)) for 
the Gaussian model. By definition we have 



R(x) = ((h(x)-h(0)) 2 ) 



Using ( p.27| ) we can write this as 



R(x) 



d 2 k e lk - 



ik-x 



i 



j J (2tt) 2 k 2 
Integrating over the angle between k and x yields 



P(k) 



R(x) = -- 



1 



d\k\ 



(J (\k\\x\)-l)P(k) 



nj Jo \k\ 

Differentiating both sides of this equation with respect to |x| yields 
1 f°° .,., 1 



(B.ll) 



(B.12) 



(B.13) 



R'(x) 



■*3 m Jo 



d\k\4(\k\\x\)P(k) = 



"./!-'' 7„ ° \\fW.U) 



where we have substituted t = \k\ \x\. For large x the integral is —1. We thus 
have 



R(x)~— In (Id) 

7TJ 



(B.15) 
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Appendix C 



Calculation of (e 



%Th) 



In this appendix we derive equation ( 3.45 ). First we derive this equation in the 
special case of a distribution T that can be written as a finite linear combination 
of Dir ac delta's. By taking appropriate limits of this special case we arrive at 

Let the distribution T be defined by: 



Th = y a n h (a 



(C.l) 



In this case we have 



(exp 



(iTh)} = exp ( iTh {1 A /exp (iTh {2 A \ 



Using (3.23) we get 

(exp (tTft< a >)) = J Dh& exp (» £„ Oi n h^ ( Xn )) 
xexp(-|/(V^ 2 )) 2 d 2 x) 

In terms of h (fc) this can be written as 



/ 



Dh™ 



exp 



E 



k L 



lev. 



( k )\ 2 ~ l ^2~/f M fc )exp 



(^/C • #nj 



(C2) 



(C.3) 



(C.4) 



Here the sum is over those k- values that belong to h^ . We write h (fc) — 
a (k) + ib (k), with a (k) and b (k) real. The integral (C.4) can be written as: 



(C.5) 



n fc ' I Ida (k) db (k) exp (-jk 2 fa {kf + b (k) 

+ 7^ J2n [«» ( a ( k ) cos ( k ■ x n) ~ b (fc) sin (fc ■ x n ))]j 

x life'/ 2 da (k) db (k) exp (- jk 2 (a (k) 2 + b (fc) 2 ) ) 

Here the prime indicates that the product is over fc- values with k y > 0. Using 



/oo 
exp(- 
-CO 



ax 



bx)dx = exp[ — )J- 



4« 



(C6) 
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the integrals over a (k) and b (k) are easily performed. (C.5) simplifies to 

2 / \ 2" 



^ «n COS (k • X n ) J + I X] "" Sin ( fc ' *") 



'(C.7) 



The product is evaluated by writing it as an exponentiated sum, and then 
replacing the sum by an integral. We thus have 



1 

Vjk' 2 



En a » C0S ( k ' X ")) + (Sn a ™ sin ( fc ' x ™)) 



rife' ex p x 

= CXp (-Y,k'vW ^n,m a n& m COS (fc • (x n - X m )) 



(C.8) 



Since the density of degrees of freedom for the field h 2 is P^ 1 (k) (see (3.15)) 
we can replace a sum by an integral as follows: 



J2-y 



d 2 k 

Jk (27T) 2 

We can thus write the exponent of ( p.8[ ) as 
1 f d 2 k 1 



P( 2 ) (fc) 



(C.9) 



2 JL2 



2j'7 (27r) 2 A: 



-f (2) ( fc ) 5Z a « a ™ cos ( fc ■ ( x « ~ *">)) 



(CIO) 



The extra factor ^ comes from taking into account the condition k y > 0. Using 



( |3.27 ) we can write this as 



~ 7 ^Ti^m^ \Xn •Em) 



(C.ll) 



from which it follows that 



(exp(lJ2n a nh(Xn))) 

exp^Y.nCtnh'-^ (x n ))exp^~^Y,n,m a nCtmG(x n ~ X m ) j ,„ { . y 



This proves ( |3.45[) in t his special case. Proving the general case is a matter of 
taking limits in ( p. 12 ). First we take the continuum limit in ( p,12[ ). Let T 
be a distribution generated by a bounded continuous function T(x) i.e. Th — 
J d 2 xT(x)h(x). Then since there exist distributions T n of the form (C.l), such 
that 



lim T n h = Th 



(C.13) 



( |3.45 ) is also valid in this case. Since for every distribution u there exist con- 
tinuous functions u n (x) such that 



lim u n h = uh 



(C.14) 



(see [|TJ), d3~45|) follows. 
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